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ABSTRACT 

New parameter free calculations including a variety of necessary kinematic and dynamic 
effects show that the results of BNL (p, 2p) measurements are consistent with the expec- 
tations of color transparency. 
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The anomalously large transmission of a hadron in the nuclear medium following or 
preceding a hard interaction is commonly referred to as color transparency [2|]. This 
topic has been actively studied theoretically ||-[(|. Color transparency depends on the 
formation of a small-sized wavepacket or point-like configuration PLC in a high momen- 
tum transfer quasi-elastic reaction. The effects of color screening or color neutrality |10 | 
suppress the interaction between the PLC and the nuclear medium. At the present ex- 
perimental energies the PLC expands as it moves through the nucleus and therefore does 
interact. This expansion is treated in Refs.|3], ITTl O. 



Color transparency is under active experimental investigation at BNL [13, [H and at 



SLAC[P1), |l6fl . The (p,2p) experiment of Carroll et al|13| finds a transparency T (ratio 
of nuclear to hydrogen cross section per nucleon after removing the effects of nucleon 
motion) with an oscillatory pattern. T increases as the beam momentum increases from 
5 to 10 GeV/c but then decreases. The SLAC experiment is unpublished at this time, but 
a preliminary report [|l(J found no effect of color transparency (CT) in (e,e'p) reactions 
for Q 2 between 1 and 7 GeV 2 . However, the stated systematic errors are about ± 10% 
which are large compared to the size of the predicted CT effects. 

The purpose of this note is to include all of the known necessary kinematic and 
dynamic effects: a proper treatment of the longitudinal (parallel to the virtual photon or 
incident proton momentum) component of the momentum of the detected protons |T7|, [18]; 
computing the measured experimental observable which involves an integration over the 
transverse momentum of the struck proton; including the effects on interference between 
PLC and non-PLC configurations produced in proton-proton elastic scattering || [7J]; and, 
a realistic treatment of the baryonic components of the PLC|l2|. Thus our aim here is 
to show that including the three kinematic and dynamic effects mentioned above along 



with the color transparency of ref. |12| leads to a natural explanation of the existing BNL 
data. No free parameters are used in the present work. 

We start the analysis by recalling the mechanisms proposed to understand the (p,pp) 
reaction. The suggestion of Ralston and Pire M , is that the energy dependence is caused 
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by an interference between a hard amplitude which produces a PLC, and a soft one (blob- 
like configuration BLC) which does not. The Ralston and Pire idea is that the BLC is due 
to the Landshoff process^]. Another mechanism with a similar effect is that of Brodsky 
and de Teramond |J in which the two-baryon system couples to charmed quarks (there 
is a small (6q) and a BLC which is a (6q,cc) object) The Brodsky-de Teramond idea is 



motivated by the fact that the mass scale of the rapid energy variation in Ann |20[ and 
in the measured transparency matches that of the charm threshold. 

Both of the mechanisms of [|7| and|| are well motivated, but imprecisely understood. 
For example, the nature of the Landshoff term is uncertain. At high energies one expects 



that this component is suppressed by a Sudakov factor plfl . At present energies, the 
energy dependence and phase of the term are not well determined, and the size of the 



object produced by the Landshoff process is not well known. See e.g. ref. p2[ . Never- 
theless, the Sudakov effects can be expected to a set of configurations with a range of 
different sizes. Similarly even if the details of the cc production amplitudes are not yet 
well- established, threshold effects will naturally lead to mixtures of BLC and PLC. Thus, 
it's natural to discuss high Q 2 elastic proton-proton scattering in terms of configurations 
of different sizes. Separating the contributing configurations into two, a PLC and a BLC 
is only a simple first step. 

The two ideas about the BLC can be used to qualitatively explain the oscillations 
observed in the (p,pp) data, but do not quantitatively reproduce features of the data Ref. 



II when combined with a proper treatment of the expansion ||11||. The agreement is 



not quantitative because including the non-zero absorption of the expanding PLC and the 
non-complete absorption of the BLC tends to make the two terms similar and weakens 
the interference effects. 

Another relevant effect is to properly account for the correct momentum of the de- 
tected protons. Consider for example, the (e,e'p) reaction in which the virtual photon 
has the momentum q. Suppose the detected proton has momentum p with p = q + k. 
For quasi-elastic kinematics p = q and k = 0. Fermi motion or final state interactions 
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can lead to a non-zero value of k. This is old news. However the non-zero components 
of k in the direction of q (k z ) have been shown to give a large numerical effect fl7l . A 



detailed study of this "Fermi motion" effect was made by Frankfurt et al [[Tj|. The mo- 
mentum k is sometimes called the momentum of the struck nucleon (and is that quantity 
if the plane wave Born approximation is valid ). We also use the term "momentum of 
the struck proton" , but only as an abbreviation. For the (p,pp) reaction the momentum 
of the detected protons can also be described in terms of the momentum of the struck 



nucleon p3 |. Indeed, the data are presented in terms of bins of k z , the component of k 
parallel to the beam direction. This is especially important for the (p,pp) reaction be- 
cause the proton- proton scattering cross section varies approximately as s -10 and here 
s = 2M 2 + 2ME LAB - 2k z P LAB « 2MP LAB (1 - k z /M) where M is the proton mass. Ref 
18| examined the (p,pp) reaction, including the effects of non-zero values of k z in their 



calculation, and obtained qualitative agreement with the 10 GeV data, but not with data 
taken at the other energies. Agreement was also obtained by Kopeliovich[24] for some 
energies. Neither group included any effects of interference between PLC and BLC. 

Our calculation of the nuclear (p,pp) cross sections uses the quantum mechanical 
treatment of color transparency developed in ref. || |^, [HJ. In that approach the time for 
the PLC to expand depends on the masses of the baryonic components of the PLC. In 
early work B |9|] only a single average state of mass M\ was used. Then the time for the 
PLC to expand is (ff^fh (M\—m) • R ecen tly we mcm ded the effects of the continuous 
baryon spectrum by using measured proton diffractive dissociation and electron deep 
inelastic scattering data to constrain the baryon masses. Our discussion of the formalism 
shall be brief and concentrate on the newer kinematic and dynamic aspects. For details 
see Refs. [g fU], |I| . 

To be definite, consider the (e,e'p) high momentum transfer process in which a photon 
of three-momentum q is absorbed and a nucleon of momentum p leaves the nucleus. As 
usual, q 2 = —Q 2 . Consider knockout from only a single shell model orbital, denoted by 
a. The observable cross sections are computed by making incoherent sums over all the 
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occupied orbitals. 

Let the amplitude be denned as Ai a , which can be described in terms of a series 
involving multiple wavepacket (PLC) -nucleon scatterings. The first term is the plane 
wave Born term B a in which the wavepacket undergoes no interactions. If full color 
transparency is obtained, this is the only term to survive. The first correction to the 
Born term is the scattering term or second term denoted by ST a . To the stated order 
we have: M. a = B a + ST a . The Born term is given by 

B a = (p\T H (Q)\a) = F(Q 2 )(p- q\a), (1) 

in which Th(Q 2 ) is the hard scattering operator. Specific effects of spin are ignored here 
and throughout this work. 

The second term is denoted by 

ST a = (p\UGT H (Q)\a), (2) 

where G is the Green's operator for the emerging small object (PLC). This object is a 
wave packet of which the nucleon is just one component. The operator U represents the 
interaction between the ejected baryon and the nuclear medium. The nuclear interaction 
can change the momentum of the ejectile and also excite or de-excite the internal degrees 
of freedom. Our approach is to treat the Green's function G as a sum of baryonic 
propagators each denoted by m. The eikonal Green's function for the emerging baryon 
is: 

9(Z - Z')e iPm( - z - z,) 
G(Z, Z>) = £ G m (Z, Z>) = J2 [ J , (3) 

m m ^Wm 

with = p 2 + M 2 — M^. Here M is the nucleon mass, M m the mass of the baryonic 
state. Our notation is that p = \p\ and p m = \pZi\- Terms beyond the first order in U 
are included by exponentiation, which is an excellent approximation for our applications 
25 . 

The work of ref. involves a simple model in which the operator U gives only one 
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excited state (m=l) when acting on a nucleon. In that model 

ST Q = -F(Q 2 ) J d 2 BdZp(B,Z) J Z ^dZ'e^ z '- z ^(Z,Z')e^ z '(B,Z'\a), (4) 

where the p(B, Z) is the nuclear density, and 

a eff (Z,Z') = a(l-e^-^ z '- z A, (5) 



and a is the proton- nucleon total cross section. This is eq (32) of ref. ||. 

Now we can exhibit the importance of the momentum (k z ) of the struck nucleon. We 
take p = q + k z for components parallel to the direction of q. Change the integration 
variable Z' to D via Z' = Z — D. Then 

/ —ik z Z , \ \ -» 

d 2 BdZap(B, Z) — ^ ^ dDe ik * D (l - e -^~^ D ) (B, Z'\a). (6) 

The effect of k z appears in two places. The e~ lkzZ factor is the same as in standard 
Born or Distorted Wave Born calculations (DWBA or Glauber optical model). The e lkzD 
factor involves a modification of the color transparency physics. To see this note that 
the real part of ST a dominates the numerics. Then cos(k z D) — cos ((k z — p + pi)D) w 
^r(p — Pi) (p — Pi — 2k z ) for small D. Since p > pi a positive value of k z reduces this 
term and increases the transparency for any given value of p. This latter effect does 
not occur in the color transparency of ref. P|, but would occur in models in which the 



baryon-nucleon interaction is treated as a finite-dimensional matrix f24], 

The above paragraph is meant as a simple explanation of a numerical effect that is 
surprisingly large. The computations of (p, pp) reactions are more involved since there is 
one incident proton wave function and two outgoing ones. Furthermore the amplitudes 
for the production of both the PLC and BLC must be taken into account. However 
the qualitative effect of including non-zero values of k z (here z is the direction of the 
beam proton) is similar to the electron scattering case. Note also that our more realistic 



calculations |0| replace the single mass Mi by an appropriate distribution of masses. 



The next step is to discuss the observables measured in ref. |T3[. Let the four- 
momentum of the target proton be denoted as (M, k). The transparency T is defined as 
a ratio T = da/da B with 
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da 

and 



f b dk z ( f dk x dk y (da A /dt)(s) (7) 



da B = J^dk z J J dk x dk y {da B / dt)(s) (8) 

in which the superscript A denotes the nuclear cross section divided by the number 
of target protons and the subscript B denotes the same quantity but computed in the 
plane wave Born approximation. In the experiment the integrations over the transverse 
momenta k x , k y were limited to about 250 MeV/c. This corresponds to almost all of the 
probability so we integrate over £tll fcjrg , ky in our calculations. The integration over the 
transverse components of k reduces the Glauber DWBA result for T by about 30%. The 
data of ref. [D| are presented in terms of T for each (k a , 



We now present the results. The calculations for an Al target (three beam momenta 
and four bins of k z ) are shown in Fig 1. The solid curves show the effect of CT using 



the complete power law form of the distribution of baryon masses as in ref. |L2[ and 
dashed curves show the results of using Ref.|| with a value of Mi = 1550MeV . This 
value gives small enough color transparency at low Q 2 so as to be consistent with the 
NE-18 data. The Ralston-Pire parameterization of da/dt for the free protons is used here 
along with their separation of the BLC and PLC terms. Details of the implementation 



of this are to be found in ref [11]. The use of the corresponding Brodsky-de Teramond 



model for the PLC-BLC interference would lead to similar results for the energy range 



we consider here |TT], pL2| . The data are from ref. fl3| . The target-dependent and target 
independent uncertainties in the normalization of T are about 10% and 25% respectively 
jl3| . We multiply the central values of the Al, Cu, and Pb data points by a factor of 0.75, 
and those of Li and C by 0.85. This is consistent with the published errors. Each data 
"point" represents a bin of k z , represented as a horizontal line. The integration of Eqs. 
(0)?® over the small bins 0.1 or 0.2 GeV/c of ref. [0 causes negligible differences with 
simply using the central values. The use of a distribution of masses starting at M + 
increases the computed da jdas at the lowest beam momentum. In either case, including 



the effects of color transparency gives good agreement with the data. 

Fig. 2, shows a summary of the different calculations for 27 Al. Here as in Figs. 1 and 
3 the solid curve represents the full color transparency calculation, the dashed curves are 
obtained by neglecting the effects of the PLC-BLC interference, and the dotted curves 
represent the use of the standard optical model or "Glauber" calculation. These latter 
curves fall far below the data, but including the effects of CT leads to a reasonable 
reproduction of the data. 

Fig. 3 shows the A-dependence of the data of ref. [|TB|. These are taken for a bin of 
k z ranging from -0.2 GeV/c to 0.1 GeV/c. Here the use of Eqs. (|7]) and (||) does matter. 
To see this compare the 27 Al data of this figure with the previous figures. The Glauber 
standard optical model leads to results (dotted curves) that again fall well below the 
data, this time for each target nucleus. The solid curves, which show our full calculation, 
are in excellent agreement with all of the data, except for the 12 GeV Al data point. The 
dashed curves show the results obtained without the BLC-PLC interference Ralston-Pire 
effect, so we see that the latter helps to account for the energy dependence, even though 
it is not a very large numerical effect. 

It is necessary to comment on the single particle nuclear shell model wave functions 
used here. Harmonic oscillator (huj = AlMeV/ A 1 ^ 3 ) wave functions are used for light 
nuclei (Li and C). The other nuclei are treated in the Hartree-Fock (HF) approximation 
with the SGII interaction of Ref. |27[| . (HF wave functions are used in obtaining the 
results of Figs. 1,2.) If the oscillator frequency is chosen appropriately, the HF wave 
functions are well approximated by a single harmonic oscillator wave function. Thus 
we find the most important effect of using the HF wavefunctions is to shift the value 
of hcu from AlMeV/A 1 / 3 to one that more precisely represents the nuclear mean square 
radius. The effect of using HF wave functions with the correct exponential dependence 
at long distance on da/da B is largest, a 20% increase independent of energy, for the 
Glauber optical model calculations for the Pb target. But the effect is negligible for the 
CT calculations. See also ref. |28| . 
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The main effect of the color transparency is to provide an increase in the predicted 
magnitude of the cross sections. One might wonder if some combination of reasonable 
effects applied to the Glauber optical model calculations yield enhancements that lead 
to reproducing the data without including color transparency effects. One possibility 



is to to claim that the experimental resolution allows the inclusion of all nuclear 
excited states in the measurement. In that case, the optical model wave functions would 
be computed using the proton-proton reaction cross section a r instead of the total cross 
section a as we have done. However, o r ja ~ 0.85 |JIJ and has only a small energy 
dependence, for the energy range of interest here. Calculations show that including this 
effect provides only 10-20 % enhancements. Including the effects of nucleon-nucleon short 
range correlations leads to a 20% increase for 12 C |28||, but is a much smaller effect for 



heavier nuclei. Thus, in the absence of CT effects, we see no possibility to elevate the 
the DWBA cross sections to the levels observed by the experiment. 

We believe the BNL experiment and the present calculations calibrate the size of color 
transparency effects. Thus we determine the energies for measurable effects. Our calcu- 
lations provide a good guide to the (e,e'p) experiments. The SLAC Ne-18 experiment is 
set for very small values of k z , so the previous results of ref |TI| stand unchanged (except 
possibly for the effects of using improved nuclear structure information for heavy targets). 
Since the predicted transparency is not large in the range of Q 2 between 1 and 7 GeV 2 , we 
expect that our earlier calculations of color transparency effects for the power-law form of 
g(Mx) will not be ruled out by the final results of the SLAC experiment fll5|, [16|| . Further- 
more, our CT calculations predict enhancements for Q 2 between 7 and 15 GeV 2 . Thus a 
higher energy experiment which examines the k z dependence should observe measurable 
effects. 

Thus, we summarize: If CT effects are included, the qualitative features of the data of 



ref [113 can be reproduced in a very natural way. No adjustment of parameters is needed. 



If the new BNL (p,pp) color transparency experiment [14] confirms the central values of 
the older results, one can be confident that color transparency has been discovered. 
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Figure Captions 



Figure 1: Effect of using distributed baryonic masses. Solid uses distributed masses. 



Dashed- M\ =1550 MeV. The data in all figures are those of Ref. [[L3[ . 

Figure 2. Full Al data. Solid- full calculation. Dashed- with out the Ralston-Pire 
effect. Dotted -Glauber. 

Figure 3. A dependence of transparency. Solid- full calculation. Dashed- with out 
the Ralston-Pire effect. Dotted -Glauber. 
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